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^ ■ Abstract 

• Three widely used scenaria for defining cuts on transverse energies of jets in ep 

. collisions are discussed. All of them are shown to suffer from the same sort of 

!>■ I unphysical behaviour when the cut regions are subject to additional constraints. 

This feature is inherent in the very way dijet cross sections are defined and cannot 
j—{ • be avoided. In particular, the symmetric cut scenario is shown to be equally suitable 

OO , for the comparison with NLO QCD calculations as the asymmetric or sum-like ones. 
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' In analyzing the data on inclusive jet production in electron-proton collisions 

>• . some lower cut Ec (at HERA typically 5 GeV) on the transverse energy 

^ , of jets is necessary for experimental as well as theoretical reasons. In parts of 

^ I the multijet phase space the cut on the transverse energies of two jets with 

- - ' highest Ex may, however, get into conflict with the infrared safety of the 

corresponding cross sections. 

In the next-to-leading QCD calculations of single and dijet cross sections there 
are at most three final state partons, which combine to two or three jets. Or- 
dering the jets according to their transverse energies in 7*p center-of-mass 
system [E^^ > E^^ > e!^^), the corresponding phase space, marked by the 
upper and lower solid straight lines in Figs. 2-4, is given by the inequality 
-S'tV^ ^ -E't'' — e!^^ . This constraint, given by pure kinematics, must then 
be combined with additional conditions that might be imposed on jets for 
experimental and theoretical reasons. Phrasing the constraint on allowed part 
of phase space in terms of E^\e^'^ does not, however, imply that only gen- 
uine dijet variables must be investigated or only the first two jets taken into 
account. For instance, the inclusive single jet cross section gets in principle 
contribution even from the third jet. 
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The aim of this letter is to argue that there is no way of choosing the cuts 
on Ef}^\ E^'^ that would be free from the sort of unphysical behaviour of cross 
sections noted in [1]. We shall illustrate our claims on the calculation of cross 
sections in ep coUisions at HERA in the region relevant for the analysis in [2] 

2 < g2 < 80 GeV, 0.1 < y < 0.85, E^ = b GeV, (1) 



but they are clearly of more general validity. All calculations reported be- 
low were obtained by means DISENT NLO Monte-Carlo program [3] using 
CTEQ6M set of parton distribution functions (PDF) of the proton. For dijet 
events there are three classes of distributions that can be measured: 



separate distributions of transverse energies and pseudorapidities 
dcr dcr dcx da 



d4')' d4'^' dr^W' d7y(2)' 



(2) 



• their sums, corresponding to distributions of what is called in [4] "trigger 
jets", 

do" d(T do" do" d(T dcr 



^Et - d^) d4'^ ' d77 - d77(i) d77(2) ' 



(3) 



• or the distributions of some of the combinations characterizing the whole 
event, hke their mean values Et = + e'^^)/2 and 77 = (77^^^ + r)^^^)/2 

da da 

(4) 



dEx ' df] ' 

but one is free to use general combinations as well. 

Symmetric cut scenario 
In this scenario the same cut is imposed on both jets 

eP > Ec, 4^ > Ec. (5) 



This is quite appropriate for the comparison of data with lowest QCD calcu- 
lations because the latter involve binary partonic hard processes which lead 
(in diparton rest frame) to two partons with the same transverse energies. Its 
use in NLO QCD calculations of dijet cross sections has been questioned in [1] 
observing that it leads to unphysical behaviour of the cross section asym{^) 
corresponding to the integral over the region defined in (5) supplemented with 
an additional constraint E^^ > Ec + A, represented graphically by the red dot- 
ted curve in the upper right plot of Fig. 1. The resulting dependence, displayed 
as red solid curve in the left part of Fig. 1, implies that for A approaching 
zero, asymiA) decreases despite the fact that the corresponding phase space 
increases. The source of this unphysical behaviour can be traced back to the 
fact that at the order aa^ perturbative QCD does not yield a well-defined re- 
sult for the double differential cross section da/dE!^^dE^^ when E^^ — E^\ 
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Fig. 1. The dependence of crfc(A) on the parameter A restricting the regions defined 
in (5,6,8) and marked in yehow, to the right of the red dotted curves. 




Fig. 2. Left: symmetric cuts scenario, Ec = 5 GeV. The soHd, dotted and dashed 
Unes correspond to fixed values of e!^^ , Et and e!^^ . As they approach Ec from 
above the lengths of the former two decrease, whereas that of Ej^' increases. Right: 
corresponding distributions of e!^^ (solid), e!^^ (dashed) and Et (dotted). 

This comes from the fact that at this order there are mass divergencies of the 
virtual corrections that must be cancelled by those of the real emissions. For 
a quantity to be calculable in perturbative QCD one must thus allow for inte- 
gration over sufficient part of the three parton phase space, where E^^ < E^j:\ 

To understand better the origin of the behaviour of asym{^) we plot in Fig. 2 
the distributions of e!^\e!^'^ and Et for E^ = 5 GeV. As expected, the 
distributions of both e!^^ and Et turn negative for e!^'^ or Et close to E^ = 
5 GeV because of the small volume (represented in left parts of Figs. 2-4 
by the length of perpendicular solid or oblique dotted lines) of phase space 
available for the real parton emission. In fact both these distributions diverge 
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Fig. 3. The same as in Fig. 2 but for asymmetric cuts scenario with Ec = 5 GeV 
and (5 = 2 GeV. 

to — oo as e!^^ or Et approach Ec from above, thereby causing agymi^) and 
o"sum(A) in Fig. 1 to bend down for small A. On the other hand da/dE^^ grows 
monotonously with decreasing A, reflecting the fact that even for E^^' = E^ 
there is plenty of phase space (represented in the left parts of Figs. 2-4 by the 
dash horizontal lines) available for the real gluon emission because eI^^ can be 
anywhere above it. Note that the integrals over all three distributions in Fig. 
2 are finite and equal as they correspond to three different ways of performing 
the integral over the same region defined in (5). There is thus no reason to 
reject the symmetric cut scenario for the comparison of data with NLO QCD, 

(2) 

but for Et close to E^ one should compare the distributions of Eq. , rather 
than those of E^^ . 



To avoid the above unphysical behaviour of agym the asymmetric cut scenario 



where 6 > 0, was suggested in [1] ^ . In Fig. 3 we plot the corresponding results 
for the same three distributions as in Fig. 2. As expected, the distribution of 
E? now behaves physically down to E^. + 6, as do those of e!^^ and E^- 
However, whereas the first distribution ends there, the other two continue 
below Ec- Both of them are, however, singular at Ej.' = Et = Ec, diverging at 
this point from left to +oo due to the divergence of real emission contributions 
as E^T^ E^t\ However, as in the case of the symmetric cuts, the integrals 
over the distributions in Fig. 3 are finite and the same. The behaviour of the 

(2) — 

distributions of Ej. and Et below Ec is of course as unphysical as the negative 
values of the distributions of Et^ and Et closely above Ec in the symmetric 
cut scenario of Fig. 2. 

^ Instead of (6) slightly different modification of (5) was used even earlier in [4]. 



Asymmetric cut scenario 



E^T >Ec + 6, Ep > E, 



(6) 
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Fig. 4. The same as in Fig. 2 but for the sum cuts scenario with Ec = 5 GeV and 
6 = 2 GeV. 

Moreover, it is simple to invent an additional constraint that cuts off part 
of the phase space region defined in (6) in a similar way as the condition 
^ > + A in the case of the symmetric cut. For instance the restriction 



A, 



(7) 



which cuts off part of the region (6) as shown in the middle right picture in 
Fig. 1 results in the integrated cross section aasymi^) exhibiting the same kind 
of unphysical behaviour for A ^ as o"sym(A). 



The cuts on e!^^ , e!^^ 



Sum cut scenario 
can also be formulated in terms of their mean value 



Et = 



(1) 



+ E?^ 



E^ > Ec. 



(8) 



As shown in Fig. 4 in this case it is the distributions in E^'^ and E^^ that are 
singular at E^. the former turns negative and diverges when E^^ approaches 
Ec from above, whereas the latter diverges to plus infinity for E^ Ec from 
below. As in the previous two cases the integrals over all three distributions 
are, however, finite and the same. 

Although even for the lowest Et the region defined in (8) involves integration 



over finite neighbourhood of the point E^^ = E^' , this scenario was claimed in 
[5] not to be infrared safe. The argument was based on the fact that similarly 
as in the case of the symmetric cut scenario (5) the integrated cross section 
asum{A), corresponding to the additional constraint e!^^ > Ec + 6 /2 + A 
behaves, as illustrated by the green dashed curve in Fig. 1, unphysically for 
A ^ 0. However, as emphasized above, the same objection can actually be 
raised against any scenario. This is clear from the geometrical meaning of the 
additional cuts involving the parameter A: as shown in the right part of Fig. 



n(2) 
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Fig. 5. The effects of the cutting the bin defined in (9-10) according to (11) on the 
distributions in e'^^ and e'>^\ 

1, in all three scenaria it cuts off the peaks of the wedges in the plane Ej:\ E^^ 
where the negative infinities of the virtual corrections dominate. 

The scenario (8) is often supplemented with a constraint on the relative differ- 
ence of jet transverse energies, for instance, {E^^^ - E^^VI^t^ + ^t^) < 1/4- 
This cuts off part of the three jet phase space but as indicated by the black 
dash-dotted curve in Fig. 1, has only small effect of the behaviour of asum{^)- 

The results presented in Figs. 1-4 demonstrate that there is no real advantage 
of the asymmetric or sum-like cut scenaria over the symmetric one. All of 
them suffer from the same kind of unphysical behaviour of two of the three 
distributions in e!^^ , e!^^ or Et and all three can be supplemented by addi- 
tional constraints that lead to the same unphysical dependence on integrated 
cross sections cr^ym(A), aasym{^), or cr^„m(A) displayed in Fig. 1. 



Distributions in bins of finite width 

In practice experimental distributions in E!p\ e!^^ and Et have finite bin 
widths. For each bin of these variables we can introduce an additional cut A, 
which chops off part of this bin, and investigate the dependence of the bin 
content on this parameter in the limit A — i> 0, much in the same way as we 
did above for the integrated cross sections. Note that for E^^ , E^^ > Ec + 6 
the choice of the cut scenario is irrelevant for these considerations as it does 
not influence the distributions there. Let us take as an example the bin 

7 < E^T^ < 10 GeV, E^T^ > E^ (9) 
7 < E^^^ < 10 GeV, E^^^ > E^ (10) 

represented in left two pictures of Figs. 5,6 by the vertical and horizontal col- 
ored bands. We first introduce the parameter A similarly as for the symmetric 
cut (6) by imposing additional constraint: 

7 + A < E^!^ < 10 GeV (11) 
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which means graphically cutting off the red areas in Fig. 5. The resulting A- 
dependence of the contributions da/dE'^'' and da/dE^^^ to the corresponding 
bin content (Tfej„(A), displayed in the right part of Fig. 5, shows that the 
distribution of e!^^ behaves unphysically for small A, whereas that of E^^ 
increases monotonously as A ^ 0. 

However, we can cut off part of the bin (9-10) in an alternative way represented 
by the red areas in Fig. 6 

7 < eP < 10 - a GeV, (12) 



which is certainly as legitimate as the cut (11). The corresponding results, 
shown in the right part of Fig. 6, lead, however, to opposite conclusion! Now 
it is the da/dErp (A) distribution that is a monotonously decreasing function 
of A, whereas da / dE!j^\A) exhibits much the same unphysical dependence on 
A as did da/dE!^\A) for the cut defined in (11). 




Fig. 6. The same as in Fig. 5 but for the cut defined in (12). 



In summary, for dijet sample there is thus no way how to introduce the cuts 
on the jet transverse energies that would avoid entirely the problem noted 
in [1]. Whatever variables, quantity and binning we choose, there is always a 
simple way how to cut the chosen kinematical region in a way that leads to 
the same sort of unphysical behaviour of cross sections shown in Fig. 1 and the 
asymmetric cut scenario is therefore not superior to the symmetric one. On 
the other hand the above exercise also suggest a simple cure: one just needs to 
choose the appropriate variable or take sufficiently wide bins in jet transverse 
energies. 
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